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ABSTRACT 

I construct a simple model of the convective burning layer during a type I X-ray burst to investigate the 
effects convection has on the stability of the layer to nonradial oscillations. A linear perturbation analysis 
demonstrates that the region is stable to nonradial oscillations when energy transport is convection-dominated, 
but it is unstable when energy transport is radiation-dominated. Thus, efficient convection always dampens 
oscillations. These results may explain the nondetection of oscillations during the peak of some X-ray bursts. 
Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries — X-rays: bursts 



1. INTRODUCTION 

Type I X-ray bursts are thermonuclear explosions that 
occur on the surfaces of accreting neutron stars. Fast 
rises and exponential-like decays lasting ~ 1 and ~ 10- 
100 seconds, re spectively, charac t erize burst lightcurves 
(for r evi ews, see iLewin et aTJll995t IStrohmayer & Bildsten 
2006). IStrohmaver et all d 19961) detected coherent oscil- 
lations in the lightcurve of a type I X-ray burst from 
the low-mass X-ray binary 4U 1728-34. Since then, as- 
tronomers have detected oscillations in burst lightcurves from 
19 additional sources (although some detections are un- 
confirmed; see lBhattacharvvall2007t iMarkwardt et al J 120071: 
IStrohmayer et al J 120081 : lLamb & B outloukosl l2008l and ref- 
erences therein). Bursts can exhibit oscillations during 
the rise, the peak, and/or the decay of their lightcurves. 
The oscillation frequency typically increases by a few 
Hz during a burst, but it asymptotes to a specific fre- 
quency unique to that source to within a few parts in 
10 3 (e.g. iGiles et aTJ 12002: Muno et alJl2002h . The asymp- 
totic frequency's stability implies that it corresponds to the 
neutr on star spin frequency (e.g. Strohmaver & Mar kwardl 
1999). The discover ies that the coherent oscillation fre- 
quency during a superburst from 4U 1636-536 matches 
the asy mptotic frequency of burst os cillations from that 
source ( Strohmaver & Markwardt 2002:), the burst oscilla- 
tion frequencies from the accretion-powered milli second pul- 
sars SAX J 1808.4-3658 (IChakrabartv et al]l200H) and XTE 



J1814-338 (Markwa rdt et all 120031: IStrohmaver et aTJ l20ljl 



match their respective known spin frequenc i es, an d the fre- 
quency of coherent pulsations ICasella et alJ (120081) detected 
in the persistent emission from Aql X-l is close to the asymp- 
totic frequency of burst oscillations from that source confirm 
this implication. 

It is most likely the rotational modulation of a growing hot 
spot on the neutro n star surface that generates oscil lations dur- 
ing the burst rise (IStrohmaver et al.llT997bl fl998l) . The time 
needed to accrete enough fuel to trigger a thermal instabil- 
ity (a few hours to days) greatly exceeds that needed to burn 
the fuel (~ seconds), which makes simultaneous ignition over 
the entire stellar s urface high l y unlikely. I gnition probably 
occurs at a point (iJossI 1 1978t ISharal [19 82). and the result- 
ing hot spot grows and engulfs the stellar surface in ~ 1 s , 
the burst rise time jFrvxell & Woosleyir i982; Bil dsten! [19951: 
IZingale et ail 120011: ISpitkovsky et alJ 120021) . 



dependent Coriolis force regulates the thermonuclear flame 
propa gation speed and t hereb y the time evolution of the hot 
spot (Spitkovsky et al. 2002; Bhattac harvva & Strohmaver! 
120071) . The accreted fuel usually ignites at the equator, and the 
resulting thermonuclear flame quickly spreads in longitude 
and generates an axisymmetric belt around the neutron star. 
No oscillations occur becau se there is no azimuthal asymme- 
try (ISpitkoyskv et alJl2002l). Ignition sometimes occurs of f 
the equator ((Cooper & Narayan 2007; Maurer & Watts 2008); 
the thermonuclear flame propagates in longitude much more 
slowly, generating a long-lived non-axisymmetric hot spot 
and hence oscillations. 

The rotationally modulated hot spot model fails to ex- 
plain oscillations during the burst decay, however, since 
the flame has engulfed the neutron star surface by this 
time, and any remnant nonaxisymmetry from ignition would 
have dissipated since the cooling time at the burst peak is 
much less than the burst duration. It is thought that ex- 
cited surfa ce modes generate nonradial oscillations in the 
burst tail ( McDermott & Taam 1987; Cumming & Bildsten 



2000t iL ee 2004; Hevl 2004, 2005; iLee & St rohmaver 120051: 
Piro & Bildstenll2005L l2006h . In particular, EeyJl 42004J) sug- 
gests buoyant r-modes as the most promising candidate, r- 
Modes travel backward in the corotating frame, so the ob- 
served oscillation frequency is less than the neutron star 
spin frequency. The r-mode frequency decreases as the sur- 
face cools during the burst decay, which explains the ob- 
served increase in the oscillation frequency (although Heyl's 
model may not adequately predict the obs erved frequency 
drift and the asymptotic freque ncy's stability; [Piro & Bildsten! 
12005b iBerkhout & Levinl 120081) . r-Modes occupy only a 
small area around the equator, so their amplitudes much be 
huge to produce the observed ~ 10% flu x variations. There- 
fore, t he modes must be driven unstable. Narava n~& Cooperl 
(2007) propose that the e-mechanism drives surface modes 
during the burst decay: if heating via nuclear burning is 
sufficiently strong relative to cooling via radiative diffusion 
and emission, some nuclear energy converts to mechani- 
cal energy and drives the oscillations (see also ILivio & Bath! 
19821: iMcDermott & Taarr] 119871; IStrohmaver & Led 119961: 
Piro & Bildsten 20041K They predict that short, powerful 
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bursts preferentially exhibit oscillations, which is generally 
in accord with observations. 

If the e-mechanism does indeed drive oscillations, then one 
would naively predict that oscillations occur preferentially 
when nuclear burning is strongest, i.e. near the lightcurve's 
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peak. This in not true; bursts exhibiting os cillations show 
them during the peak only about half the time dGalloway et alJ 
2008). Some bursts are so powerful that the nuclear lu- 
minosity temporarily exceeds the local Eddington luminos- 
ity; the extra thermal energy converts to both kinetic energy 
and gravitational potential energy and expands the outermost 
layers of the neutron star. Astronomers often detect oscil- 
lations from these so-called photosph eric radius expansion 
(PRE) bursts dMuno et alj|2000ll200lh . bu t they rarely detect 
oscillations during the PRE phase itself (jSmith et al] 1 19971; 
Strohmaver et alJ Il997at iMuno et all 120021; iGallowav et all 



2008). Why, then, do bursts often fail to exhibit oscillations 
during the peak phase? 

To address this question, one must consider the mechanisms 
by which the burning layer transports energy. Highly effi- 
cient convection is the primary energy transport mechanism 
during the firs t few seconds of t he initial phase of a power- 
ful burst (e.g lJosslll977L [T978] lHanawa & Sugimoto)[l98l 
Wallace et all 1 19821 iPaczynskil Il983t IWoosley et al l 120041 
Fisker et aLll2008l) . It is well known that convection can either 



drive or dampen pulsations in variable stars (e.g.[Coxlll980t 
IStellingweril 19861: lUnno et alj|1989t iDupret et al.ll2005l and 
references therein). Does conv e ction affect burst oscillations 
as well? iMcDermott & Taaml d!987l) were the first to sug- 
gest such a coupling between convection and pulsations dur- 
ing type I X-ray bursts, but they did not pursue this idea. In 
this investigation, I develop an analytical one-zone model of 
the convective burning region during a type I X-ray burst to 
determine the effects convection has on the pulsational stabil- 
ity of the burning layer. I find that efficient convection damp- 
ens oscillations, much like it does in stars near the red edge of 
the instability strip. These results may explain the nondetec- 
tion of oscillations near the peaks of some type I X-ray bursts. 
I begin in ^2] by deriving the equations that govern both the 
oscillatory and convective motions. I perturb the governing 
equations and derive the stability criterion in $3] an d I a Pply 
this criterion to type I X-ray bursts in §|4| I conclude in §0 
with a discussion of the results. 

2. OSCILLATION AND CONVECTION EQUATIONS 

Star ting from the work of lUnnol (1 19671) and iGabriel et alJ 
( 1 1975b . I construct a one-zone model of the convective burn- 
ing layer during a type I X-ray burst on a rapidly rotating 
neutron star. My goal is to conduct a linear stability anal- 
ysis on the entropy equat ion in a manner similar to that in 
iNaravan & Cooper! d2007l) and determine the stability of the 
configuration to nonradial oscillations. Usually, one conducts 
such calculations assuming the unperturbed configuration is 
static. This assumption is inapplicable to the problem at hand, 
for even the unperturbed convection zone involves vertical 
motions of the matter. Thus it is necessary to distinguish be- 
tween fluid motions due to convection and fluid motions due 
to oscillations. A further complication is that an entirely sat- 
isfactory theory of convection does not yet exist. I use the 
standard mixing length theory to describe convective motions. 
Specifically, I demand that the mixing length theory of con- 
vection applies to both the initial configuration and the per- 
turbation. Constru cting such a model is not trivial (see, e.g. 
Gabri eTet al.lll975l) . Luckily, the parameters relevant to type 
I X-ray bursts allow me to simplify the problem enormously. 

I make the following assumptions: 
(I) The Boussinesq approximation of the mixing length theory 
applies to the conve ctive motions (see, e.g. lCox & Giulill968t 
lHansen et al.ll2.Q04b . In particular, I ignore pressure fluctua- 



tions in the entropy equation and density fluctuations in the 
continuity equation. 

(II) The mixing length I is much smaller than the typical non- 
radial oscillation wavelength. This assumption is easily justi- 
fied, for the mixing length is presumed smaller than the ver- 
tical pressure scale height (~ meters), whereas the oscillation 
wavelength is of order the neutron star radius (~ 10 km). 

(III) The lifetime of a convective element r c is the short- 
est relevant timescale other than the vertical sound crossing 
timescale. Specifically, r c <C f2 _1 , ui~ l , r nuc , and tr, where Cl 
is the neutron star spin frequency, u> is the mode angular fre- 
quency, T nuc is the nuclear burning timescale of matter within 
a convective element, and tr is the radiative cooling timescale 
of a convective element. The last inequality implies that con- 
ve ction is efficient. Es t imate s of r c range from 10~ 6 -10~ 3 
s (IMcDermott & Taaml 119871: ICumming & Bildstenl 120001: 
IWeinberg et al.ll2006l: iFisker et al.ll2008l) . in accord with this 
assumption. 

(IV) Both the turbulent pressure and the rate of turbulent ki- 
netic energy dissipation into heat are negligible. 

(V) The layer is always in vertical hydrostatic equilibrium. 

I proceed as follows. First, I write down the hydrodynamic 
equations. These equations describe the aggregate behavior of 
both the oscillatory and convective motions. The convective 
motions have short wavelengths of order the mixing length /, 
which is much smaller than the oscillation wavelengths by as- 
sumption (II). The disparity of the two lengthscales allows me 
to separate the equations governing the oscillatory motions 
from those governing the convective motions. In this work, 
I investigate only the stability of the oscillations to nonadia- 
batic perturbations, not the properties of the o scillations them- 
selve s (for a discussion of the latter, see. e.g. | Longuet-H iggins 
1968: iBildsten et alJfl996l: INaravan & Coopedl2007l) . There- 
fore, I conduct a linear stability analysis on only the entropy 
equation governing the oscillatory motion. The perturbed en- 
tropy equation contains perturbations of quantities describing 
the convective motions. I solve for the perturbed convective 
quantities in terms of quantities describing the oscillatory mo- 
tions by perturbing the convective equations. 

2.1. Fundamental Hydrodynamic Equations 

The fundamental continuity, momentum, and entropy equa- 
tions of the model are 



dp 
dt 



+ V-( P v) = 0, 



(1) 



d\ 1 

— + (v-V)v + 2fixv=-ra--Vn, (2) 

at p 



(3) 



where p is the density, v is the velocity, g is the gravitational 
acceleration, p is the pressure, T is the temperature, s is the 
entropy per unit mass, e is the nuclear energy generation rate, 
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-vr 4 



(4) 



is the radiative flux, a is the radiation constant, and k is the 
radiative opacity, g is presumed constant since the vertical 
thickness of the convective layer is small relative to the stel- 
lar radius. Equations (HO govern the aggregate behavior of 
the burning layer during a type I X-ray bursts. Specifically, 
they describe both the oscillatory and convective motions of a 
matter element. To proceed further I need to derive both the 
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equations that govern only the long wavelength, oscillatory 
motions and the equations that govern only the short wave- 
length, convective motions. I do so in the following subsec- 
tion. 

2.2. Oscillatory and Convective Equations 

Mot ivated by the work of lUnnol (1 19671) and iGabriel et all 
(1 19751) . I decompose the physical quantities into their oscilla- 
tory and convective parts: 

p=p+Ap, 
p=p+Ap, 
T = T + AT, 

s = s+ As, 

v = u + V, (5) 

where p, p, f, s, and u are values averaged over a horizon- 
tal area with dimensions much larger than the mixing length 
but much smaller than the oscillation wavelength, and Ap, 
Ap, AT, As, and V are the local quantities that describe the 
convective fluctuations. The convective quantities are much 
smaller than their respective oscillatory quantities except for 
the velocity, in which case |V| 3> |u|. In fact, u = in the un- 
perturbed configuration. Note that Ax = for all quantities x. 
Furthermore, I define the Lagrangian derivative following the 
oscillatory motion 



d d 
dt dt 



(6) 



For each of the three fundamental equations (Q][3]l, I first de- 
rive the oscillatory equations by taking the horizontal average 
of the corresponding fundamental equation. I then derive each 
convective equation by subtracting the corresponding oscilla- 
tory equation from the fundamental equation. 

2.2.1. Continuity Equation 

Implementing equations (f5]l and ©, I write the continuity 
equation (fl]i as 

d(p+Ap) 



dt 



■ + ( ( o+Ap)V-u + V-(pV) = 0. 



(7) 



Setting Ap = by assumption (I) and notin g that pV = ( since 
convection involves no bulk motion, e.g. lGabrielHl996l) . the 
space-averaged continuity equation becomes 

dp 



dt 



+ pV-u = 0. 



(8) 



Subtracting equation © from (Q~|i and again implementing as- 
sumption (I) (Ap = 0, p = constant) gives the continuity equa- 
tion describing the convective motions, 

V-V = 0. (9) 

2.2.2. Momentum Equation 

Taking the horizontal average of equation (f2]l gives 



du 



1, 



, +(V-V)u+(V-V)V+2f2xu = --Vp-gz. (10) 
dt p 

Subtracting equation ( TTOb from (f2]) gives the momentum equa- 
tion for the convective motions, 

dV 



— + (V • V)V- (V ■ V) V+ (V • V)u - (V • V)u 

dt 

+2Vtx\ = --V(Ap) + ^-Vp. 

p p 2 



(11) 



To be consistent with the mixing length theory, I set 
(V-V)V-(V-V)V=A— , 



(V- V)u-(V- V)u = A— , 



where 



\V;\ 



(12) 



(13) 



is the lifetime of a conve ctive element, V 7 = V ■ z, and A is 
a number of order uni ty ( Bohm- Vitensd fl958t lUnnol [19671 : 
Griga hcine et alJfeOOSl) . The momentum equation for con- 
vection becomes 

dV 1 Ap u V 

— +2ftxV= — V(A/?) + — ^Vp-A A-. (14) 

dt p p L T c T c 

2.2.3. Entropy Equation 
Taking the horizontal average of equation (0 gives 



-ds 



pT— = pe-V -F R -pTY ■ Vs. 
dt 



(15) 



I need to express the last term in a more useful form. Using 
the thermodynamic identity Tds = dH-dp / p, where H is the 
enthalpy per unit mass, I find 



pTV • V* = pV-VH-V- Vp. 



(16) 



V • Vp is the turbulent kinetic energy dissipation rate and is 
negligible by assumption (IV), so I omit it from the subse- 
quent analysis. Assumption (II) implies that VH = V(H + 
AH) w S7AH. Using this and equation ©, I find 



pTY ■ Vs = V ■ pVAH. 



(17) 



In general, AH = T As + Ap / p by the thermodynamic iden- 
tity, but since Ap = by assumption (I), AH = TAs, so the 
horizontally-averaged entropy equation becomes 



pf^-=pe-V-(¥ R + ¥ c ), 
dt 



F c = pVrAs, 



(18) 
(19) 



where Fc is the horizontally-averaged convective flux. I use 
equations ( TT8l and ( fT9] l in $3] Subtracting equation ( fT5l ) from 
(0) convective entropy equation 



ds 



dAs 



( ApT + pAT) — + pT —— + pTX ■ Vs - pTX ■ Vs 
dt dt 



= pe-pe-V-(F R -F R ). (20) 

To simplify the equation and keep it consistent with the mix- 
ing length theory, I se t, in analogy with equations (TT2l see 
Grigahcen eet al.il2005l) . 



V • A(pT Vs) - V • (pT Vs) = pT 



As 



(21) 



Setting A(pe) = and AFr = by assumption (III) gives the 
convective entropy equation 



Ap AT\ ds dAs wr „ As 
—^- + ^) —+ + V- Vs = . 

p T J dt dt t c 



(22) 
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3. PERTURBATION OF THE HEIGHT-INTEGRATED EQUATIONS 

For simplicity, I construct a one-zon e model for the verti- 
cal structure of the accreted layer (e.g. [ Fujimoto et al.llT981h 
lBildsten|[l998t iNaravan & Cooped 12001 ^ I assume that the 
horizontally-averaged quantities describing the accreted mat- 
ter are constant throughout the layer. By assumption (V), I 
write the height of the accreted layer as 

P 



The equation of state is 

pk^T aT 4 

P= + -5- 

pm p 3 



pm p 



In 



PS 



+ constant, 



(23) 



(24) 



(25) 



where k B is Boltzmann's constant, /i is the mean molecular 
weight, m p is the proton mass, and (3 = pk B T j pm p p is the 
ratio of gas to total pressure. Perturbing the equation of state 
gives 

ST _ 1 Sp _ 1 Sp 



XT P V T P 



(26) 



where xt = 4 - 3/3 and 



( dlnp 



4-3/3 



(27) 



\dlnTjp 

is minus the coefficient of thermal expansion. From equations 
©, ( fl8l l, and ([T9i l, the height-integrated entropy equation is 



_-ds 
pT di" 



■pe-g(F R +F c ), 
acgT 4 



(28) 



(29) 
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F c = pTV z As. (30) 

The goal of this work is to perturb equation d28l > and deter- 
mine the stability of the burning layer to nonradial oscilla- 
tions. However, the expression for the perturbed convective 
flux SFq includes perturbations of convective quantities. To 
proceed, I first need to express the perturbed convective quan- 
tities SV Z and S As in terms of the perturbed oscillatory quan- 
tities Sp, Sf, and Sp. The following subsection is devoted to 
this. 

3.1. Perturbed Convective Equations 

I conduct a standard WKB perturbation analysis on the con- 
vective equations. I assume that all dynamical quantities vary 
as exp(-iu!t + /k ■ x), where k = k x X + k y y + k z i and kh <C 1 . 
Since the lifetime of a convective element r c is assumed to be 
much smaller than the mode period 2tt/ui, I presume that all 
unperturbed convective quantities are in quasi-steady state. 

Deriving the perturbed continuity equation for convection 
is simple. From equation (0, the continuity equation for con- 
vective motions is 

k-V = 0, (31) 

so the perturbed continuity equation for convection is 

k-<5V = 0. (32) 

Invoking the WKB approximation, integrating over z, and 
simplifying, the momentum equation for convection ( fl4l i be- 
comes 

^ + 20xV = -^ + ^z-A^-aI, (33) 

dt p p T c T c 



where I have assumed that V/> ~ {dp/dz)i in the unperturbed 
state. Perturbing equation d33l , setting both d/dt — > and 
£1t c — » by assumption (III), and ignoring the term Au/r c 
because | V| ^ |u| gives the perturbed convective momentum 
equation 



Pg J gTc\ T C J V P 



(34) 



Equation ( |34l contains the perturbed convective quantities 
SAp, 5Ap, <5V, and St c , whereas the desired perturbed quan- 
tities are 5 As and 5V Z . To remedy this, I proceed as follows. 
Taking the dot product of equation d34i l and k and using equa- 
tions (ED and {32]!, I find 



(35) 



Taking the dot product of equation ( 1341 ) and z and using equa- 
tion d35T l then gives 

f(t-t)-mH). - 

Conducting a similar procedure on the unperturbed, steady- 
state convective momentum equation, I find 



It follows that 



gTc \P J\ k* 



SV 1= 5t 1 S(Ap/p) 
V z r c + Ap/p ' 



(37) 



(38) 



To proceed further, I need expressions for the two terms on the 
right hand side of the above equation. From equation dT3b . 



St c _Sl _ SV Z 



(39) 



I make the standard assumption that the mixing length / is 
proportional to the pressur e scale height h (e.g., ICox & Giulil 
119681; lHansen et al] |2004|) . The lifetime of a convective 
element is negligible compared to the oscillation period 
(i.e. u>t c -C 1) by assumption (III). Therefore, I assume the 
mixing length instantaneously adjusts to changes in the scale 
height. It follows that 



(40) 



SI Sh Sp Sp 
I h p p ' 

where the last equality follows from equation ( l23l . Note that 
equation d40b would not apply if, for example, ujt c 3> 1; in 
that case, SI /I = 0. Equations d39b and ( f40b then give 



St c Sp Sp SV Z 

Tc P P V z 



(41) 



Writing the convective entropy as TAs = c p [AT - 
(dT /dp) s Ap], setting Ap = by assumption (I), and 
noting that AT/T = -(1 jvj) Ap/p from the equation of state, 
I get 

— =-v T — , (42) 
P c P 



where 



Cn = 



5k B Z32-24/3-3/3 2 



p 2pm p 



5(3 2 



(43) 
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5 



is the specific heat at constant pressure. Thus, I find that 

5V, 1 / Sp Sp SAs Sv T Sc„ 



V, 



1 / Sp Sp SAs Sv T 

2 \ p p As vt 



— ■ (44) 



All I need now is an expression for SAs in terms of the oscil- 
latory variables. 
Integrating the convective entropy equation d22i > over z 



gives 



An AT \ ds dAs V z spg As 

"^ + Tb; + ir + ^ = — ■ (45) 

p 1 J dt dt p t c 



Perturbing this equation and setting wr c 
(III), I find 



by assumption 



5V Z Ss Sp Sp SAs 5t c 
V z s p p As t c 



(46) 



Using equations (|4TT > and (04), the above relation simplifies to 

SAs _ 5s 
As ~ s ' 



(47) 



3.2. Perturbed Entropy Equation 

Perturbing the left side of the horizontally-averaged entropy 
equation ( TT~5T > and using equation d26l i to eliminate ST/T, I 
find 



5 pT 



Pgh 
(T 3 -l) 



S_p 
P 



■r, 



(48) 



where Tj = (dlnp/dlnp) s and I^ — 1 = (dlnT /dln p) s are the 
usual adiabatic exponents (e.g. lHansen et al.ll2004l) . The nu- 
clear energy generation rate e is a function of both density and 
temperature. For small perturbations about the initial config- 
uration, I write 



dine \ ( <91ne 

d\nT J ' \dlnp 7 T 
(49) 

where ftj is the characteristic heating time of the layer via 



. / ST Sp 
Se=il h gho v-=r + r]— 
T p 



nuclear burning. Using equation ( l26l i to replace ST/T in favor 
of Sp/ p, the perturbed heating term becomes 



5e = n h gh 



v Sp ( v \ Sp 
-+ V — 

Xt p \ v T 



(50) 



The layer cools via both radiative diffusion and convection. 
The radiative flux Fr depends on the radiative opacity k, 
which is a function of both density and temperature. For small 
perturbations about the initial configuration, I write 



5 k 



Kip 



dlriK 
d\nT 



<91nK 
d\np 



(51) 

Note that equations (148115 1|> are simply generalizat ions of their 
analogous equations of lNarayan & Coop er j2007). The mag- 
nitude of the unperturbed convective flux Fq depends on the 
convective quantities V, and As. I cannot determine their val- 
ues in terms of the oscillatory quantities p, f, and p from first 
principles, so the magnitude of Fc is uncertain. I parameterize 
this uncertainty by defining the dimensionless quantity 



fc 



Fc 



Fr + Fq 



(52) 



the ratio of the unperturbed convective flux to the unperturbed 
total flux. Using equations $29$, CP, (ED, and (|52j, I find 



5(F R + F c ) = n c phUl-f c ) 



(4-C)f-?-^ 
T p p 



'Sp ST SV Z SAs 
p T V z As 



(53) 



where is the characteristic cooling time of the layer. Fi- 
nally, using equations ( l44b and d47l ) to eliminate the perturbed 
convective terms and using equation d26l ) to eliminate ST/T, 
the perturbed cooling term becomes 



5(F R + F c ) = n c ph 

T4-C 

Vt 



4-C 

XT 



3 3-C 



2 vt 



Xt 



Sp 
P 



Sp 

p 



2 \ s vt 



Sep 
cp 



(54) 

Substituting equations d48b . d50b . and d54l > into the perturbed 
entropy equation 



■ S(pe)-gS(F R +F c ) (55) 
and simplifying gives and expression for Sp in terms of Sp, 

(56) 



5p =A Sp 



P P 

where A is a dimensionless complex quantity. For purely adi- 
abatic perturbations, A = Ty . In the limit where the convective 
layer is gas pressure dominated (J3 — » 1) 

5 / 1 + (2i/5)(fic/w)[4 - C+ e/c(7/2+ £- - $l h /n c )(v-r,)] 



A = - 



1 • ' i+(2i/3)(n c /w)[3-c-/ c (3/2-o-(n«/n c )(i+i/)] 

(57) 

and where the layer is radiation pressure dominated (f3 — > 0) 

i + a/4)(n c /ioM+ M l /2-Q+ (n h /n c )r,] 

3 V 1 + (i/3)(n c /uj)[-(/4 + / c (3/2 + C/4) - (n h /n c )( l + 

(58) 

I work exclusively in these two limits hereafter for simplicity. 

4. PULSATIONAL STABILITY 

I invoke the quasi-adiabatic approximation to determine the 
linear stability of the layer. That is, I presume the fractional 
entropy change during one oscillation period is small. An 
equivalent statement is that the heating and cooling rates are 
small relative to the oscillation frequency, i.e. f2/,/w <C 1 
and il c /uj <C 1. The typical heating and cooling times dur- 
ing a type I X-ray burst are of order one second, so SI/, 
and f2 c ~ Is" 1 , where as the mode s have angular frequencies 
~ lOOradianss" 1 (e.g., Hevl 2004), so the quasi-adiabatic ap- 
proximation is valid. 

I use the following criterion to determine the pulsational 
stability of the convective layer: it is unstable if the work done 
on the layer during one oscillation period is positive, and it is 
stable if the wor k done on the layer during on e oscillation pe- 
riod is negative (Cox 1980t IDnno et al.|[l989t) . The pressure 
and density perturbations both vary in time as exp(-/ajf), but 
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they are in general out of phase because A is complex. If 
Im(A) > 0, Sp lags behind Sp and the area traced in the P— 

V diagram during one oscillation period is positive; the PdV 
work done on the layer is positive and the layer is unstable. 
Conversely, if Im(A) < 0, Sp leads Sp and the area traced in 
the P—V diagram is negative; the PdV work done on the layer 
is negative and the layer is stable. For purely adiabatic per- 
turbations, A is real and Sp and Sp are in phase. The time 
evolution of Sp and Sp traces a curve of zero area in the P— 

V diagram, so the PdV work done in one cycle is zero; the 
modes neither grow nor decay. Thus the sign of Im(A) deter- 
mines the stability of the layer. 

At the high temperatures achieved during a burst, the opac- 
ity scales accord ing to the approximate formulae given by 
iPaczvnskil (119831). which g ive -0.1 < £ < 0, -0.5 < C < 0. I 
follow iNara van & Cooper! (120071) and set £ = and C = -0.25 
for simplicity. Solving for the imaginary part of A in equa- 
tions (l57l i and ( l58l . I find the instability criterion for nonradial 
oscillations to be 



7 + 5/ c 



n c 10+677+4^ 

when the layer is gas pressure dominated (/3 — > 1), and 



> 



1 + 17/e 



n c 16+ 12?7+4zy 



(59) 



(60) 



when the layer is radiation pressure dominated ([3 — > 0). Note 
th at equatio ns ( 1591 and (l60l reduce to equations (76) and (77) 
of iNaravan & Cooper! (120071) when f c = 0, as they should. If 
the heating rate due to thermonuclear burning during a burst 
is sufficiently large relative to the cooling rate, some of the 
nuclear energy converts to mechanical energy and drives the 
surface modes. It is clear from equations 09l and ( f60T > that 
convection dampens oscillations in both cases. For 0=1, 
the effect is fairly modest; the minimum heating rate needed 
for instability when cooling is convection-dominated is less 
than a factor of 2 greater than that needed when cooling is 
radiation-dominated. However, the temperature of the layer 
during a burst is often > 10 9 K, especially when a sizabl e 
convective layer develops dJossI 119771: IWoosley et al.l 12004). 
so the case j3 = is more relevant for this work. For j3 = 0, the 
effect is huge; the minimum heating rate needed for instabil- 
ity when cooling is convection-dominated is 18 times greater 
than that needed when cooling is radiation-dominated. There- 
fore, oscillations during type I X-ray bursts are unlikely when 
c onvection dominates the cooling. 

iLed (120041) draws the opposite conclusion: nuclear burning 
drives oscillations in a convective layer. Focusing primarily 
on the mode structure instead of the stability^ criterion, Lee 
simplifies his calculation by assuming <5(V • Fc) = 0, which 
implies that the convection zone is "frozen" during the pertur- 
bation. This assumption is inapplicable to type I X-ray bursts 
because lot c <C 1, i.e. the convection zone has sufficient time 
to respond to oscillatory motions. 

5. DISCUSSION 

In this investigation, I have found that efficient convection 
dampens nonradial oscillations during type I X-ray bursts. 
This may explain the nondetection of oscillations near the 
peak of some bursts. The basic physics of this effect is simple. 



Consider the lateral compression of a column of matter during 
a burst. When the matter is compressed, the temperature and 
density rise and generally increase the nuclear energy genera- 
tion rate; the cooling rate may either increase or decrease. If 
the marginal increase in the heating rate is sufficiently large 
relative to the cooling rate at maximum compression, the pres- 
sure will contin ue to rise and thereby dri ve the oscillations 
(see also §5.2 of Naravan & Cooper 2007). Energy transport 
by convection is more effective than that by radiative diffusion 
when the layer is compressed, so the cooling rate is larger at 
maximum compression if the layer is convective. In this case, 
the pressure is more likely to decrease at maximum compres- 
sion and hence dampen the oscillations. 

The convection zone lasts only a few seconds after ignition 
and disappears roughly when the burst luminosity reaches its 
peak. The e-mechanism may drive oscillations at this time. 
The r-mode growth ti mescale is a few seconds wh en the burn- 
ing layer is radiative (Narayan & Cooper 2007), so oscilla- 
tions may be detected in the burst tail even if convection 
dampens them near the peak. 

It is unclear if convection alone explains the paucity of 
burst oscillation detections during PRE. Strong convection 
is guaranteed during a PRE burst since the nuclear burning 
timescale is much less than the radiative cooling timescale, 
but the heating rate is probably large enough relative to the 
convective cooling rate to satisfy the instability criterion. It 
may be that, although convection cannot damp the oscilla- 
tions entirely, it lowers the growth rate enough to keep the os- 
cillation amplitudes below the detectability threshold. Alter- 
natively, perhaps the surface modes are unstable during PRE, 
but the bursting layer's vertical extent is sufficiently large rela- 
tive to the mode wavelength to smear out the nonaxisymmetry 
dCumming & Bildstenll2000l) . 

I have presented a simplified, analytical, one-zone model of 
the convective burning layer as a first attempt to study the ef- 
fect convection has on the driving of burst oscillations. There 
are several important issues I am unable to address. First, 
convection alters both the energy transport mechanism and 
the temperature gradient. I have investigated only the for- 
mer in this work. Second, the model assumes the convection 
zone's extent is constant in time. During an actual burst, the 
convection zone first grows to lo wer pressur e s and then re- 
cedes to the ignition region (e.g., Joss 1978; Woosle v"et al.1 
120041 IWeinberg et al.ll2006l iFisker et al]|200l . Presumably, 
the convection zone in a given column of matter will expand 
when the column is compressed and thereby further dampen 
the oscillations, but this is merely a speculation. Third, I 
cannot determine the actual values of the critical parameters 
Oft/fi c , / c , v, and 77 of the stability criterion. One can deter- 
mine th em only with detailed, mul ti-zone calcu l ations such as 
those of IWoosley etail d2004b and lFisker et all d2008l) . 
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